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Infinite, neutron-proton symmetric, neutral nuclear matter has a critical temperature of Tc = 
17.9±0.4MeV, a critical density of pc = 0.06±0.01 A/fm^ and a critical pressure of pc = 0.31±0.07 
MeV/fm'^. These values have been obtained from our analysis of data from six diflerent reactions 
studied in three different experiments: two "compound nuclear" reactions: ^**Ni+^^C— ^-^''Se and 
^''Ni+^^C— ^'^''Se (both performed at the LBNL 88" Cyclotron) and four "multifragmentation" re- 
actions: 1 GeV/c TT+^^''Au (performed by the ISiS collaboration), 1 AGeV ^"^Au-h^^C, 1 AGeV 
"'^La-h^^C and 1 AGeV ^''Kr-h^^C (all performed by the EOS collaboration). The charge yields 
of all reactions as a function of excitation energy were fit with a version of Fisher's droplet model 
modified to account for the dual components of the fluid (i.e. protons and neutrons). Coulomb 
effects, finite size effects and angular momentum arising from the nuclear collisions. 
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I. INTRODUCTION 

In the long history of the study of the liquid to vapor 
phase transition of nuclear matter [Illl4j various studies 
have sought to determine one or more critical exponents 
[U m fTOHT^ [Ti] . other studies have examined caloric 
curves [S], and others have reported the observation of 
negative heat capacities [9j. All of these efforts suffer 
from the lack of knowledge of the systems location in 
pressure-density-temperature (p, p, T) space. Specifi- 
cally, interpretations of caloric curves and negative heat 
capacities depend on assumptions of either constant pres- 
sure or constant density (TSl US] • In the case of determin- 
ing critical exponents, it was assumed that the systems 
were at coexistence and that the surface energy was the 
single, dominant factor. The analysis presented below 
makes no assumptions about the location of the system 
in (p, p, T) space and accounts specifically for other ener- 
getic considerations (e.g. the Coulomb force and angular 
momentum) . 

Our approach begins with the time honored idea of 
considering nuclei as drops of a hypothetical nuclear fluid. 
The liquid drop model 17] takes up this idea quantita- 
tively. The approximately constant binding energy per 
particle in heavier nuclei suggests that this fluid is bound 
together by a saturating short range force similar to that 
acting between the molecules of simple fluids (i.e. Van 
der Waals like). 

Present day formulations of the liquid drop model 
[IHl US] express the binding energy in terms of a volume 
term proportional to the number of nucleons A and cor- 
rected for finiteness by means of an expansion in terms 
of A^^i'^ of which only the first (surface energy) order 
term is kept. Additional corrections are added to ac- 
count for neutron/proton asymmetry. Coulomb interac- 
tions and pairing effects. 

Global fits to nuclear masses lead, on the one hand 
to a reproduction of binding energies to within 1% (<10 
MeV), and on the other to the characterization of the 
hypothetical fluid mentioned above, where finiteness. 



neutron/proton asymmetry and Coulomb have been re- 
moved. This is the "bulk nuclear matter" which has been 
studied theoretically over the history of nuclear physics 

m- 

Van der Waals fluids admit various phases, among 
which are the vapor and the liquid. Thermodynamically, 
the equilibrium coexistence of these phases and the asso- 
ciated liquid-vapor phase transition are well understood. 
The Van der Waals aspects of the nuclear binding energy 
lead naturally to the question: does nuclear matter sus- 
tain a vapor phase as well as the condensed liquid phase? 
Is there in the phase diagram a (first order) coexistence 
line terminating at a critical point jT]? If so, how can 
one obtain such information experimentally? 

We answer those questions (in the affirmative) and de- 
termine the coexistence curve and critical point of bulk 
nuclear matter from data obtained in three experiments 
and six different reactions. 

In this paper, we describe first a physical picture of 
the nuclear reactions in question. Then we provide a 
brief description of the experiments. Greater detail about 
the experiments can be found in the references provided. 
Next, we give a detailed description of the theory and 
analysis used on the experimental data. Finally, we use 
the results of that analysis to determine the coexistence 
curve and critical point. 



II. THE PHYSICAL PICTURE 

Thermal nuclear sources (compound nuclei and higher 
energy nuclear aggregates) emit particles such as neu- 
trons, protons and heavier charged fragments into vac- 
uum in a process that is very similar to evaporation |21| . 
This type emission from thermal, equilibrated systems 
is in contrast to the direct, or prompt, particle emission 
from excited nuclear systems out of equilibrium. 

For fluid systems like water, evaporation rates allow 
one to recover the properties of the saturated vapor in 
equilibrium f22|. However, in nuclear systems, finiteness 



and the presence of Coulomb effects prevent such a simple 
approach to the characterization of the phase diagram. 

It has been shown previously how it is possible to "re- 
duce" the nuclear evaporation rate to that of an infinite, 
uncharged symmetric fluid: finiteness is accounted for in 
terms of the "complement" approach [23], Coulomb ef- 
fects can similarly be factored out "M] , and the corrected 
rates can then be related to the properties of the hypo- 
thetical nuclear matter vapor [22) . 

We now demonstrate explicitly how properties of the 
bulk nuclear matter such as its phase diagram can be as- 
certained from experimental measurements of fragment 
distributions starting with a physical picture of frag- 
ment production from excited nuclei. This illustrates 
how one can talk about coexistence without the vapor 
being present and how an equilibrium description, such 
as Fisher's theory [23 [5S] (described in section IV A), 
is relevant to the free vacuum decay of an evaporating, 
nuclear system [53]. 

Thermodynamicians would determine a phase diagram 
via direct measurements of the pressure, density and 
temperature of their fluid. However, such direct mea- 
surements of temperature, density and pressure for a nu- 
clear fluid are problematic. On the other hand, the mea- 
surement of clusters in nuclear reactions has been easily 
achieved and has a long tradition. Thus, we believe that 
in nuclear physics, this is the royal avenue toward the 
extraction of the phase diagram. 



A. The virtual vapor 

Let us now consider a liquid in equilibrium with its 
saturated vapor. At equilibrium, any particle evaporated 
by the liquid is restored on the average by the vapor 
bombarding it. In other words, the outward evaporation 
flux from the liquid to the vapor is matched by the inward 
condensation flux. This is true for any kind of evaporated 
particle. The vapor acts as a mirror that reflects the 
evaporated particles back into the liquid. 

One could probe the saturated vapor by putting a de- 
tector in contact with it. However, since the outward 
and inward fluxes are identically the same at equilib- 
rium putting, the detector in contact with the liquid also 
probes the vapor. Therefore, we do not need the vapor 
to be physically present in order to characterize it com- 
pletely. We can study the evaporation of the liquid and 
dispense with the surrounding saturated vapor. In these 
situations one thinks of a virtual vapor, realizing that 
first order phase transitions depend exclusively upon the 
intrinsic properties of the two phases, and not on their 
interaction. Of course, if the vapor is not there to restore 
the emitting system with its return flux, evaporation will 
proceed and the system will cool. 

An excited nucleus is a small drop of equilibrated nu- 
clear matter that emits neutrons, protons and higher 
charged fragments into vacuum according to statistical 
decay rate theory. In this picture there is no surround- 



ing vapor, no confining box and no need for either. As 
described in the preceding paragraph, by studying the 
outward flux of the first fragments emitted from a ther- 
mal source at equilibrium, we can study the nature of the 
vapor even when it is absent (the virtual vapor). 

Quantitatively, the concentration nji^j{T) of any species 
A{ constituents at temperature T in the vapor is related 
to the corresponding decay rate RAt{T) (or to the decay 
width r^f(T)) from the nucleus by matching the evapo- 
ration and condensation fluxes 

RaAT) = ^-^Y^ « flAAT) {vAAT)MnAvA,)) , (1) 

where VAf{T) is the thermal velocity of the species A{ (of 
order yjT j Ai) crossing the nuclear interface represented 

by the cross section Uinv (of order ^s where A^ is the 
mass number of evaporating nucleus), T is the temper- 
ature of the equilibrated, excited nucleus when the first 
fragment is emitted, h is Planck's constant and Coulomb 
effects have been, for the time being, neglected (they are 
dealt with below). 

Equation (fl]) shows the fundamental and simple con- 
nection between the (compound nucleus) decay rate and 
the fragment concentration in the vapor. Thus, the va- 
por phase in equilibrium can be completely characterized 
in terms of the decay rate. 

The physical picture described above is valid instanta- 
neously. The result of successive evaporation in vacuum 
leads to abundances of various species of emitted frag- 
ments that arise from a continuum of systems at different 
temperatures J27| . This leads to complications in various 
thermometers: kinetic energy, isotope ratios, etc. 

Our way of avoiding this complication is to consider 
only fragments that are emitted very rarely so that, if 
they are not emitted first, they are effectively not emitted 
at all. In other words, we consider only fragments that 
by virtue of their large surface energy (and high charge) , 
have a high emission barrier. The rapidly increasing 
Coulomb barrier with fragment charge Z strongly en- 
hances this effect. Thus a lower cut-off of about Z « 6 
is used in the analysis that follows. 



III. EXPERIMENTS 

The above physical picture is used here to analyze the 
data from two kinds of experiments: compound nuclear 
decay and multifragmcntation. Both types of experi- 
ments measure the total yield or number of fragments 
emitted from a thermal nuclear source, ^^^(-E*), as a 
function of the excitation energy of the source, E^ , and 
charge of the fragment, Zf. For both types of experi- 
ments it is assumed that the collisions produce an ex- 
cited, equilibrated thermal source of radius rg consisting 
of As nucleons [Z^ protons and N^ neutrons) at excita- 
tion energy E* and with angular momentum /. This is 
the initial state of the system: an excited, thermal nu- 



cleus which emits neutrons, protons and heavier charged 
fragments. 



A. Compound nucleus 

The first kind of experiment gives rise to a compound 
nucleus [SSj. A compound nucleus is formed when one 
nucleus impacts another nucleus and the two combine to 
form a single, compound system. The nucleon number 
and charge of the compound nucleus is just the sum of the 
nucleon number and charge of the two colliding nucleii. 
Its excitation energy can be determined from the energy 
of the bombarding nuclei and the masses of the target and 
projectile. The excited compound nucleus is a thermal 
source that emits protons, neutrons and other heavier 
charged fragments. 

The compound nucleus experiments analyzed here 
were performed at the 88" Cyclotron of the Lawrence 
Berkeley National Laboratory ^29]. An Advanced 
Electron-Cyclotron-Resonance (AECR) ion source [30] 
was utilized to produce highly charged ^^Ni and ^"^Ni ions 
which, after injection into the cyclotron and acceleration 
to the desired energy, impinged on a high purity [31] car- 
bon target (1.0 mg/cm^). The fragments emitted in the 
reactions were detected in two position-sensitive AE — E 
detector assemblies placed on either side of the beam. 
The methods of the energy and position calibrations of 
the /S.E and E detectors have been described previously 



B. Nuclear multifragmentation 

The second kind of experiment analyzed in this work 
gives rise to a phenomenon called multifragmentation 
[71 I33H39] . In a multifragmentation experiment, one nu- 
cleus is accelerated to a high velocity and impacts an- 
other nucleus and in the experiments considered here, 
one of the colliding nuclei is larger than the other. Typ- 
ically, the collision between nuclei in multifragmentation 
experiments is more violent than the that in compound 
nucleus experiments. The two nuclei either partially fuse 
or a "fireball" is generated from the occluded parts of 
the target and projectile. The larger of the two nuclei 
promptly loses nucleons during the collision leading to an 
excited, thermal remnant with a smaller nucleon number 
and charge than the initial nuclei. In the experiments 
considered here, the greater the excitation energy of the 
remnant, the the smaller the nucleon number and charge 
of the remnant. Figure [l] shows the nucleon number, 
As, and charge, Zg, as a function of excitation energy, 
E* , for the remnants created in the multifragmentation 
experiments considered here. The excited remnant is a 
thermal source that emits protons, neutrons and other 
heavier charged fragments. In multifragmentation exper- 
iments, the excitation energy is estimated based on mea- 
surements of the kinetic energy of the fragments emitted 
from the remnant and other considerations |40l [41] . 
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FIG. 1: Top: the nucleon number as a function of excitation 
energy for the thermal sources created by the multifragmenta- 
tion experiments. Bottom: the charge of the thermal source 
as a function of excitation energy. 



Some of the multifragmentation data analyzed here 
are from the reactions 1 AGeV ^^'^Au-I-^^C, 1 AGeV 
"^La-fi^C and 1 AGeV ^'^Kr+'^'^C and were collected 
by the EOS Collaboration at the Lawrence Berkeley Na- 
tional Laboratory Bevalac. This experiment studied the 
projectile fragmentation and detected nearly all of the 
charged reaction products on an event-by-event basis 
[H |40H42] . Charged particles with charges of from 1 to 
6 were identified using a time projection chamber [43] 
while a multiple sampling ionization chamber detected 
charged particles with charges from 7 to 79 [44] . 



2. IStS 

The other multifragmentation data analyzed here are 
from the reaction 1 GeV/c tt-I-^^^Au and were collected 
by the ISiS Collaboration at the Alternating Gradient 
Synchrotron (AGS) at Brookhaven National Laboratory 
[^, 45-48.. The AGS provided beams of 1 GeV/c tt 
incident on a gold target. Particles with charges from 1 to 
16 were measured by the Indiana Silicon Sphere (ISiS) 47r 
detector array [45| providing a high statistics, exclusive 
data set with finer bins in excitation energy than the EOS 
experiment [46 . 



IV. ANALYSIS 

We now provide a derivation of a formula for the aver- 
age fragment yields that should be observed in the exper- 
iments based on the initial state of the system (the ex- 
cited, compound nucleus or remnant) and the final state 
of the system (the fragment and its complement). We 
start from Fisher's droplet model [HI [5B] and modify it 
to account for effects that arise from finite size [23|, the 
Coulomb force |24j . isospin, angular momentum and the 
secondary decay of excited fragments. 



A. Fisher's droplet model and the complement 

Fisher's droplet model [551 [55] is a physical cluster the- 
ory that has successfully: described the cluster distribu- 
tions in percolating systems [42] and lattice gas (Ising) 
systems j49] ; reproduced the compressibility factor at the 
critical point |50j : predicted (within a few percent) the 
compressibility factor of real fluids from the triple point 
to the critical temperature [SH [52] ; and has been used to 
describe the nucleation rate of real fluids [S31 [S3] . 

Physical cluster theories of non-ideal fluids assume that 
the monomer-monomer interaction is exhausted in the 
formation of clusters, and that the resulting clusters be- 
have ideally (i.e. they do no interact with each other) 
[55|, I56j . Further, clusters of a given number of con- 
stituents A{ can be characterized by a chemical potential 
(per constituent) fi and a partition function qAf(T,V) 
that depends on the temperature T and volume V of the 
fluid and is given by 



qAf{T,V) 



/ 2nmAfT 



exp 



AG 



(2) 



where V is the volume, niAf is the mass of a fragment of 
A{ constituents. Here AG is the free energy cost for the 
formation of that cluster [57l 



AG^AE-TAS + pAV 



(3) 



where AE and AS are the energy and entropy cost of the 
formation of the cluster respectively and p is the pressure 
and Ay is the change in volume due to the formation of 
the cluster. 

Because of the ideality of the fluid of clusters, the pres- 
sure and density are readily determined. The pressure p 



T 



P= V E'?^.(^'^)^ 



At 



V 



(4) 



At = l 



and the density p is 

-. CO 



(5) 



where z is the fugacity z = e'^/'^ . The concentration of 
size Af clusters is then 



nAt{T,z) 



qAAT,V)z^^ 
V 

^At f2TrmA/rY 

\ /l2 



go exp 



AG 

T~ 

AfAfi\ f AG\ 
^^)exp(- — j 



(6) 



where A/i is a measure of the distance from coexistence 
in terms of the chemical potential which, following Fisher 
[251 , absorbs the thermal wavelength and qo is a normal- 
ization constant. At coexistence, A^ — 0, the cluster 
concentration is given by 



UAt (T) = qo exp 



AG 



(7) 



There have been many derivations of AG, but here we 
follow a general derivation using the complement method 
[23] and concentrate on the change in free energy between 
the initial state and the final state. Because in the com- 
plement derivation bulk terms do not survive, they are 
omitted in the derivation below. 

The initial state consists of an equilibrated liquid drop 
consisting of As particles. In the final state the drop 
has just emitted a cluster or droplet or fragment with A{ 
particles. Also in the final state along with the fragment 
is the complement. The complement is what is left of 
the drop after the fragment has been emitted and thus 
consists of ^c = ^s — Af particles. 

In determining the free energy of the initial and final 
state, we follow Fisher's contribution to physical cluster 
theory which was to endow clusters with a surface energy 
and to provide an estimate for the entropic part of the 
free energy associated with the formation of a cluster 
[251126]. 

Since the vapor of clusters is ideal, its internal energy 
is given by 



A=A, 



E= J2 "^aEa 



(8) 



A=l 



where Ea is the binding energy of a cluster and is deter- 
mined via a liquid-drop expansion 



Ea 



a„A + a, A" 



(9) 



where a^ is the bulk of volume energy coefficient, a^ is the 
surface energy coefficient and a is an exponent describing 
the relationship between the surface and volume of the 
fragment. One intuitively expects that a ~ 2/3 for three 
dimensional systems. 

The only contributions to the change in the energy 
between the final and initial states are those that are 
associated with the change in the net surface area. Thus 
the energy change is given by 



At = l 



AE- 



[{As - AiY + A\ 



AZ 



(10) 



Similarly, the entropy of the fluid is given by 



S = '^ uaSa- 



(11) 



A=l 



Fisher conjectured that the entropy of a cluster could 
also be estimated with a liquid drop type expansion 



Sa = byA + bsA'' -tIuA 



(12) 



where &„ is the bulk of volume entropy coefhcient and bs 
is the surface entropy coefficient. See reference jJS] for 
more on the origins of the logarithmic term. 

As with the energy, the only contributions to the 
change in the entropy between the final and initial states 
are those associated with the change in the net surface 
area. Thus the entropic contribution is given by 



AS = b^iA ~ AfY + A^ ~ a:] 
-iA,-At)A{- 



r In 



A. 



(13) 



The contribution to the free energy due to the change 
in volume, the pAV term is negligible compared to the 
energetic and entropic contributions and will be ignored 
here. 



Equations (pi) and (111 combine to show that the 
change in free energy is 



AG = {a,^Tbs)[{A,-Afr + A^^A-J 
■{As-At)At- 



+ Trln 



^s 



(14) 



At the critical temperature, Tc, the surface's contribu- 
tion to AG vanishes leaving only the logarithmic term, 
thus eq. ( 14 1 shows that Tc is defined as 



Tr.= 



(15) 



Assuming that the coefficients a^ and bs are indepen- 
dent of the temperature the quantity (a^ — Tb^) can be 
rewritten as a,£ with 



Now the fragment concentration becomes 

-{A,-Ai)Ai~ 



(16) 



UAfiT) = go 

X exp 



As 

a^e [{As 



A,r + A1^Al] 



T 



In the limit of a bulk liquid where Ag — > oo, 
As — v4f w As and we obtain 



UAi (T) = qoAf ^ exp 






(17) 
then 

(18) 



which is the expression Fisher derived for the cluster con- 
centration [25 1 126 1 157 ] 

At the critical point, e — (the surface free energy 
vanishes) and the fragment distribution is given by a 
power law. The power law has been explicitly verified 
in percolation and Ising systems |6l [lOl IH |Ml I59H8T] 
and implicitly verified in a wide variety of physical fluids 

[51 IS]. 

For the present work, the liquid is not inflnite and at 
most As « 175, so equation (18) cannot be used and all 



the terms in equation (17) must be used, in addition to 



other terms that arise due to the nuclear nature of the 
fluid, e.g. a Coulomb term, an isospin term, etc. Those 
other terms are derived below by examining the initial 
and final states of the evaporating, equilibrated, excited 
nuclear source. 



B. Characterization of the initial state: Properties 
of the thermal source 

The nucleon number of the excited, equilibrated, evap- 
orating nuclear source is Ag with Zs protons and Ng 
neutrons. For the compound nucleus experiments, the 
source is defined as the sum of the target and projectile. 
For the multifragmentation experiments the source was 
measured in the experiments and found to be a nucleus 
smaller in nucleon number than the larger of the projec- 
tile or target as shown in Figure ITj 

The temperature of the thermal source T is estimated 
via the Fermi gas. The excitation energy of the source in 
terms of MeV per nucleon is related to the temperature 
of the source by 



E: = aT^ 



(19) 



where the level density parameter, a, is modified to ac- 
count for the empirically observed change with excitation 
energy [52] and is given by [53] 



1 



1 



( AE*s 

l^£;bind 



(20) 



In Eq. ([20| E^'"'^ is the binding energy of the ther- 
mal source. Because the fragment yield distributions 
analyzed below are measured as a function of fragment 
charge {Z{) and because the fragment mass {A[) is only 
estimated, pairing and shell effects are neglected and the 
binding energy of a nucleus of nucleon number A {Z pro- 
tons and N neutrons) is found (in MeV) from the liquid 
drop expansion [T7HI9] 



rpbind 
J^A,Z 



-ay (l - hy'^) A 
a, (1 - hy^) A" 

3Z{Z~l) 
"^5 roAi/3 



(21) 



where a„ 
1.0884, K 



15.7335, K = 1.6949, a^ = 17.8048, k^ = 
1.43997 MeV fm, tq = 1.2181 fm and the 



asymmetry (or isospin) parameter is 

A-2Z 

y^ — -A — ■ 



(22) 



The values of the parameters are taken from reference 

m. 



We note that Eq. (21) already gives us as, one of the 



parameters needed to determine the critical temperature 
as shown in Eq. (15 1. 



In the case of the compound nucleus experiments, the 
angular momentum I of the thermal source was esti- 
mated |29j . However, we multiplied that estimate by 
a constant, /q, that was left as a fitting parameter. For 
the multifragmentation experiments, we parametrized / 
as 



i/o+/i^: 



(23) 



where the coefficients /q and /i are also left as fit param- 
eters 

The energy of the thermal source due to its angular 
momentum is (classically) 



El = 



T-sr^ 



(24) 



All radii in this work are taken to be 

r = roAi/3. (25) 

Following Fisher, the entropy of a nucleus is estimated 
I with a and r set to their three di- 
ass value: a = 0.63946 ± 0.0008 and 



based on Eq. (12 



mensional Ising c 

T = 2.209 ±0.006 gg. 

The free energy of the thermal source is the free energy 
of the initial state 

Ginitial = Eg™ + Eg —TSs- (26) 

The pV contribution to the free energy is neglected here. 

C. Final state: the fragment and complement 

The final state considered here is immediately after the 
emission (or evaporation) of a neutron, proton or heavier 
charged fragment. The fragment has mass mi, radius ri 
and Af nucleons [Zf protons and N{ neutrons). 

After fragment emission the thermal source is reduced 
in nucleon number and labeled as the "complement" and 
has mass rric, radius Tc and A^ nucleons (Zc protons and 
Nc neutrons). It is assumed that the the fragment and 
complement are both spherical, at normal nuclear density 
and that the surface of the fragment and the surface of 
the complement are in contact. 

Conservation of angular momentum dictates that the 
fragment and complement system have the same angular 



momentum as the thermal source. The energy associated 
with this angular momentum is then (classically) 



ELc = -. ?■ (27) 

'+^ |(^^^2+„,^^2) + 2^^J^(re + rf)^ 

The Coulomb energy between the fragment and com- 
plement is 



E, 



Coulomb 

f+c 



ZsZc 



Tf-t-Tc 

The free energy of the final state is 



(28) 



G 



final 



= E, 



bind 



E. 



bind 



E, 



f-l-c 



E, 



Coulomb 
f-l-c 



- T(5f-t-S'c) 



(29) 



The determination of quantities associated with the frag- 
ment and complement are discussed below. 



1. Properties oj the fragment 

The charge of the fragment Zf is measured by the de- 
tectors in the experiments. The nucleon number of the 
observed fragment A^^^^ is estimated via the EPAX pa- 
rameterization 84 



^EPAX 



-0.10167 + 1.9638Zf + 0.0057221Zf2. (30) 



The excitation energy of the fragment is estimated 
with the Fermi gas relation 



E'i = aAfT^ = a (A 



EPAX 



iVcvap) T' 



(31) 



where Ai is the nucleon number of the fragment prior 
to any secondary decay or evaporation (under the ap- 
proximation that only neutrons are evaporated from the 
fragments). 

The number of neutrons A'ovap that can be evaporated 
from a nucleus {^A^^^"^ ,Zi) is approximately: 



iV,, 



E* 



B„ + 2T 



(32) 



where -B„ is the neutron binding energy of the nucleus in 
question. This is estimated as 



B„ 



E 



bind 

^EPAX_i 2 



- m 
E 



.Zt 



bind 



Combining Equations (31) and (32) gives 



N,, 



^j,2^EPAX 



B„ 



2T - aT^ 



(33) 



(34) 



where ^2 is a fit parameter to account for the crude na- 
ture of this approximation. 



The fragment's initial nucleon number A{ is 



The inverse cross section for fragment emission is 



Af^A 



EPAX 



N, 



cvap- 



(35) 



This estimate assumes that only neutrons are emitted 
during the secondary decay. 

The fragment's binding energy (Ef^'^'^), radius (rf) and 
entropy {S[) are determined using equations (21) and 
(12), respectively, with A{ and Zf. 



2. Properties of the complement 

Conservation of mass and charge give the mass ^c, 
charge Z^ and neutron Nc number of the complement. 
The binding energy {E^™'^), radius (re) and entropy (Sc) 
of the complement are determined using equations (21) 
and (12 1, respectively, with A^ and Zc- 



D. Average fragment yields 



(Tinv = 47r (rf + rcY 



(39) 



where only the geometric cross section is considered since 
the Coulomb effects are explicitly dealt with below [24) . 
Equation (ItI) shows that the concentration of the vir- 
tual vapor depends on the free energy cost of cluster for- 
mation AG. For a fragment emitted from an excited 
nucleus AG is given by equations ( 26 ) and ( 29 1 and is 



AG 



ZTibind 



E. 



bind 



^f+c 



TT'Coulonib 



E' 



bind 



- Tn4A?+A^,-A^)-T\n 



AfAc 
A. 



-Ei 
(40) 



This can be simplified and written as: 



AG — Gfinal — G^ 



final 

asAI - 
AGoo 



litial 



T{b,A", 



T In Af ) -I- A/infs 



(41) 



All of the experiments discussed here measure the aver- 
age yield of fragments with a given charge as a function of 
excitation energy of the thermal source Yz^ (E*) where: 



yzAe:) = 



number of Z[ fragments in events with E* 
total number of events with E* 



and its associated error on the mean 6Yzi {E* ) . 

It is assumed that there is a one to one relationship be- 
tween the excitation energy of the thermal source and the 
temperature and a one to one relationship between the 
charge of a fragment and its nucleon number. Thus the 
fragment charge yields as a function of excitation energy 
are equivalent to the fragment mass yields as a function 
of temperature which will be written as YAf,Zi (T). 

For the first fragments emitted their yield is given by: 



Ya.^z, (T) = AtRA,,z, (T) = Y 



(36) 



where At is the time duration of the measurement of the 
decay of the thermal source and Rai.Zi (T) is the rate of 
fragment emission. 

Equation (II]) shows that the fragment emission rate 
from the thermal source is related to the concentration 
'iT'Ai.Zt (T) of any species (^f , Z{, ) in the "virtual" vapor 
which matches the evaporation (or emission) flux out of 
the thermal source with the condensation flux into the 
thermal source [22] 



RAs.Z! (T) « (wAfCinv) nAf,Z! (T) . 



(37) 



The mean thermal velocity of the fragment normal to 
the plane of emission is given as 



[VA!) = 



T 



2iTmAi 
where an ideal vapor has been assumed. 



(38) 



where A/infs absorbs all the terms in eq. (40) (and is 
divided by Af) not explicitly written in eq. (41 1. A/ijjfs 



is an "effective chemical potential" that arises due to the 
nuclear nature and finite size of the thermal source. The 
factors not absorbed in A/Xnts describe the free energy 
cost of the formation of a fragment from bulk nuclear 
matter which is written as AGoo. 



E. Fitting the experimental charge yields 



TABLE I: Fit details 



Reaction 


Points 


Number 


of 


Zt 


Et (AMeV) 




fit 


Parameters Range 


Range 


''»Ni+^"C^'"Se 


54 


3 




[6, 16] 


[1.13,2.02] 


^*Ni-fi2C^''''Se 


40 


3 




[7,151 


[1.08,1.82] 


1 AGeV ^-^Kr-h^'C 


26 


4 




[6, 13] 


[1.75,4.75] 


1 AGeV "^La-h^^C 


53 


4 




[6, 181 


[1.75,4.75] 


1 AGeV is"Au-h"C 


96 


4 




[6,25] 


[1.75,4.75] 


1 GeV/c TT+^^^Au 


234 


4 




[6, 15] 


[1.50,4.00] 



The fragment charge yields are given by 

YAt,Zt (T) = At (ufCTinv) go exp f - -—- j 



At (ufCTinv) exp 



A/infs 



T 



X (Jo A exp ( ^ 



(42) 



a formula that depends on several quantities that are 
unknown, e.g. the At, q^, etc. However, the ratio of 



TABLE II: Fitting results 



Reaction 


xi 


d.2 


/o 


h 


bs 


Tc (MeV) 


b«Ni+^^C^™S6 


1.3 


0.1 ±0.1 


1.20 ±0.09 


- 


0.97 ±0.02 


18.4 ±0.3 


64Ni+12c^76sg 


0.4 


0.5 ±0.2 


1.3 ±0.2 


- 


0.99 ±0.01 


18.0 ±0.2 


1 AGeV ***Kr+''C 


3.3 


0.0±5 X 10~ 


-5 -83 ±13 


18 ±4 


1.02 ±0.01 


17.5 ±0.2 


1 AGeV "^La+^^C 


1.1 


1.8 ±0.1 


19 ±15 


-18 ±3 


0.973 ±0.008 


18.3 ±0.2 


1 AGeV ^^Uu+^'C 


1.3 


1.1±0.1 


3 ±23 


7±6 


1.007 ±0.007 


17.7±0.1 


1 GeV/c TT+^^'^Au 


3.2 


0.0±3x 10- 


"' 151.1 ±0.7 


-3.8 ±0.2 


1.032 ±0.001 


17.26 ±0.02 
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FIG. 2; The critical temperature as a function of the thermal 
source mass. The results of the compound nuclear reactions 
are shown with circles and the results of the multifragmen- 
tation reactions are shown with squares. Colors show the 
results for different experiments. The solid line shows the av- 
erage of all the measurements and the dotted lines show the 
RMS variation: 17.9 ± 0.4 MeV. 



the yield of a fragment of a given charge at a given ex- 
citation energy YAf.Zf (T) to some reference yield of a 
fragment with another charge at the same excitation en- 
ergy Yji^'z'. {T) cancels the all constants of proportional- 
ity and several unknown quantities. Therefore, the ratio 
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FIG. 3: The angular momentum values for all six reactions 
predicted by the fitting of the charge yields. Error bars are 
smaller than the points for the compound nucleus and ISiS 
experiments. The dashed line shows a thermal estimate of 
the angular momentum [85]. Dotted lines show an estimate 
of the angular momentum based on the colliding nuclei in the 
multifragmentation experiments. See text for details. 



of charge yields was fitted with the reference yield taken 
as the yield of fragments with the charge equal to the 
lower limit of the Z{ fit range listed in Table |T] The ratio 
of the charge yields is given by 



Ya,,z, (T) _ 

YA'^,zi (T) {v'.aU 

/AG 
X exp 



exp ' '"'' 



V T 



T 



Af ^ exp 



aAfe 
T 



e^rexpi-^ 



(43) 



where Q absorbs all factors other than those in Fisher's 
model for the bulk fluid. The only free parameters in 



eq. (43 1 are ^2, the secondary decay coefficient, Iq (and /i 
for the multifragmentation data) the angular momentum 
parameter (s) and b^, the surface entropy coefficient. For 



most of the fragments considered here, secondary decay 
results in the evaporation of < 3 neutrons. 

Previous attempts similar to the fitting of data sug- 
gested by Eq. (431 have been made. Reference [3] fit in- 



clusive fragment yields with a form of Fisher's model that 
used the liquid-drop model to parameterize the energy 
cost of fragment formation. However, that work did not 
account for finite size effects, the fragment-complement 
Coulomb energy or the effects of angular momentum. In 
references [TTl|31j the multifragmentation data presented 
here was fit with a version of Fisher's model that included 
an ad hoc Coulomb energy parameterization, but did not 
account for the effects of finite size or angular momen- 
tum. In reference |22) some of the compound nucleus 
reaction data presented here was fit with Fisher's model 
where all the effects of finite size. Coulomb energy and 
angular momentum were absorbed into a chemical po- 
tential which yielded much less physical insight than the 
present analysis. 

Table IT] shows the fragment charge and excitation en- 
ergy range over which the fits were performed. The 
lower limit in the fit range of Z^ for the reactions 
^''Ni-|-^^C->^^Se was set by the available data while for 
the other data sets it was for fragments sufficiently large 
to ensure that they were emitted first, or not at all. 
The upper limit in the fit range of Zi for the reactions 
58]^j_l_i2(j_j.70gg g^j^j TT-l-i^'^Au is determined by the avail- 
able data while for the other data sets it is determined by 
the largest fragment for which the fragment, complement 
scheme is appropriate, i.e. Zi < Z^jl or the largest value 
of Zi present in the data, whichever is smaller. The range 
in excitation energy is determined by starting at excita- 
tion energy values where shell effects cease and where 
there are one or fewer fragments in the Zi fit range. 

For the compound nucleus data, there are three fit pa- 
rameters for each data set: 6, di and /q. For the multi- 
fragmentation data, there are four fit parameters for each 
data set: 6, ^2, -^o and I\. On average, there are ~ 23 
points per fit parameter. 

The results for &, /q, I\ and ^2 are fisted in Table [iTl 
Figure |2] shows that the results for the critical tempera- 
ture of bulk nuclear matter determined from these exper- 
iments as a function of the mass of the thermal source. 
The multiple points shown in Fig.[2]for the EOS and ISiS 
experiments are due to the differing mass and charge of 
the excited thermal sources; the masses and charges of 
the excited thermal sources are shown in Figure [ll 

Figure [3] shows the result for the angular momentum 
from the fits. Also shown in that figure are other esti- 
mates of the angular momentum. One estimate is ther- 
mal in nature and is based on the angular momentum 
that would be imparted by the evaporation of a source 
beginning with 7 = |85| . The thermal estimate serves 
as a lower limit on the estimate of the angular momen- 
tum. The other estimate is based on the impulse and 
impact parameter of the collision; this serves as an up- 
per limit on the estimate of the angular momentum. The 
results for the estimate of the angular momentum from 
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FIG. 4; The fragment yields from the ^^Ni-f^^C-^^Se data. 
The curves show the fit to the data. There are 54 data points 
that are fit with three free parameters. 



EOS multifragmentation data lies, to within error bars, 
between the two estimates. 

The results of the surface entropy coefficient hs from 
all the experiments agree to within 3%. Combining this 
estimate of hg with the value of a^ (given below Eq. (21 )) 



gives an estimate of the critical temperature of bulk nu- 



clear matter as T^ = 17.9 ± 0.4 MeV via Eq. (15). This 
value agrees well with theoretical predictions |371 



Figure [2] shows the value of T^ as a function of the mass 
and charge of the thermal source. 

Figures [i] through [9] show the fragment charge yield 
ratios as a function of the fragment charge. Figures [l0| 
through |15| show Arrehnius plots in the form of the frag- 
ment charge yield ratios as a function of the inverse tem- 
perature. In all of those figures, the data is shown by 
the empty circles with error bars and the fits are shown 
with solid lines. The solid lines are segments drawn (to 
guide the eye) between the fit values at each Zi or 1/T. 

Alternatively, one can combine all the results shown 
in figure g through |T5| by plotting Ya,,z, {T) IYa',^z[ (T) 
divided by QA^^ as a function of aA^e/T. This col- 
lapses the all measured fragment yield ratios for any Af, 
Zi and E* onto a single curve. This is shown for all 
the data sets in Figure [16] Because the finite size ef- 
fects. Coulomb effects, etc. have been scaled away (via 
dividing the yield ratios by 8), these plots show the co- 
existence curve of bulk nuclear matter in terms of the 
concentration of droplets of the bulk nuclear fiuid that 
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FIG. 5: The fragment yields from the '^"'Ni+^^Q^^'^Sc data. 
The curves show the fit to the data. There are 40 data points 
that are fit with three free parameters. 



FIG. 7: The fragment yields from the 1 AGeV ^^^La+^^Q 
EOS data. The curves show the fit to the data. There are 
53 data points fit with four free parameters. 
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FIG. 6: The fragment yields from the 1 AGeV '^^Kt+^^C EOS 
data. The curves show the fit to the data. There are 26 data 
points fit with four fit parameters. 



FIG. 8: The fragment yields from the 1 AGeV ^''''Au+^^C 
EOS data. The curves show the fit to the data. There are 
96 data points fit with four free parameters. 
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FIG. 9; The fragment yields from the 1 GeV/c tt+^'^'^Au ISiS 
data. The curves show the fit to the data. There are 234 
data points fit with four free parameters. 



FIG. 11: The Arrehnius plots from the ''"'Ni+^^C-^^Se data. 
The curves show the fit to the data. There are 40 data points 
fit with three free parameters. 
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FIG. 10: The Arrehnius plots from the ^^^Ni+^^C^^Se data. 
The curves show the fit to the data. There are 54 data points 
fit with three free parameters. 



FIG. 12: The Arrehnius plots from the 1 AGeV '^''Kr+^^C 
data. The curves show the fit to the data. There are 26 data 
points fit with four free parameters. 
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FIG. 13: The Arrehnius plots from the 1 AGeV ^^^La+^^C 
data. The curves show the fit to the data. There are 53 data 
points fit with four free parameters. 



FIG. 15: The Arrehnius plots from the 1 GeV/c vr+^'^^Au 
data. The curves show the fit to the data. There are 234 
data points fit with four free parameters. 
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FIG. 14: The Arrehnius plots from the 1 AGeV ^^AU+^^C 
data. The curves show the fit to the data. There are 96 data 
points fit with four free parameters. 



would comprise a saturated vapor in equilibrium with an 
infinite, bulk nuclear liquid. 

In the plots shown in Figures [4] through 16 all of the 
data (all fragments of all charges and all excitation ener- 
gies) for a given experiment were fit simultaneously. The 
circles show the data points (errors are shown when they 
are larger than the symbols). The solid curves show the 
fit to the data. 



V. CONSTRUCTING THE PHASE DIAGRAM 

The fitting of the data as illustrated above gave the 
critical temperature of bulk nuclear matter. Once T^ is 
determined, it is possible to determine the entire coex- 
istence curve of bulk nuclear matter which completely 
maps the liquid-vapor phase diagram. 

The first step is to determine the coexistence curve in 
reduced units: 



— , — and — . 



Pc Pc 



Tn 



(44) 



where p is the pressure, p is the density and the subscript 

"c" denotes values of the quantities at the critical point. 

It is assumed that the formation of fragments exhausts 

all non-idealities, so that the pressure and density can be 
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FIG. 16: The scaled charge yields for all six reactions. 



obtained by simple sums. The pressure is 
p = Tj2nA{T) 

A 



Tj2<loA-^e^p 



T 



and at the critical point 



A A 



The density is given by 

= ^go^^^''exp[- 

A ^ 

and at the critical point 



T 



A A 



(45) 



(46) 



(47) 



Using the reduced quantities removes the unknown nor- 
malization qq. All other quantities in the above sums 
are known. The errors associated with Tc^ t and a are 
propagated to generate errors on the reduced quantities. 



A. Reduced density 

The empty squares in Figure [17] shows the vapor 
branch of the p-T phase diagram of nuclear matter, al- 
beit in reduced form. Those points were constructed by 
performing the sums in equations ( 47 ) and ( 48 1 . 



The empty circles in Figure [T7| show the liquid branch 
which was determined as follows. First, Guggenheim's 
universal function describing the reduced pi.v/Pc-T/Tc 
phase diagram [SD] 



Pl.v 



1 + die ± dpe^ 



(49) 



(48) (where e is given by eq. (16)) was fit to the empty squares 



on Figure 17 from 0.55Tc < T < Tc which is roughly the 
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FIG. 17: The reduced density-reduced temperature coexis- 
tence curve for bulk nuclear matter. Empty squares show the 
vapor branch. Empty circles show the liquid branch. Solid 
curves show the results of the fit to the vapor branch. Dot- 
ted curves show the extrapolation of that fit. The dashed line 
shows the extrapolation of the law of rectilinear diameter. See 
text for details. 



range over which Guggenheim's function describes dozens 
of fluids: from the triple point to the critical point. Here 
di and dp are left as fit parameters and /? is a critical 
exponent and is [551 IE] 



/3 = 



= 0.3265 ±0.0001. 



(50) 



The vapor branch is described by the eq. ( 49 1 with the 
minus sign and the liquid branch is described by the 



eq. (49 1 with the plus sign. 



The solid curve on Figure [17| shows the result when the 
empty squares (vapor branch) were fit with eq. ( 49 1 which 
resulted in di = 0.04315 ± 0.00001 and dp = 1.15714 ± 
0.00001; the errors quoted are those resulting from the 
fitting procedure. These values are different from those 
that Guggenheim found, but that is true for other fluids 
as well, e.g. helium and mercury [52] • However, eq. (49) 



still describes the coexistence curve of those fluids, albeit 
with different values for di and dp. 

For the fluids that Guggenheim examined, the liquid 
branch (for the range 0.55Tc < T < Tc) was described 
by equation ( 49 ) with the sign of dp changed. The solid 
curve shown with T/T^ > 1 on FigurellTlshows the result. 



Dashed curves on Figure 17 show extrapolations for 
T < Q.55Tc. The extrapolation for the vapor branch 
shows unphysical behavior with Pv/pc < for T/Tc < 
0.25, thus some care must be taken when determining the 



pi/ Pc-T/Tc coexistence curve at low temperatures. The 
Pv/ Pc-T/Tc coexistence curve at low temperatures has 
already been determined from equations (|47|) and (|48 1 . 



To determine the liquid branch of the coexistence curve 
for low temperatures we start with the the law of recti- 
linear diameter |90 which is 



Pi + Pv 
2pc 



1 + die 



(51) 



We extrapolated this linear function in e from T ~ 0.55rc 
to T = 0. This is shown by the dashed line in Fig- 
ure [ITJ We then used that extrapolation and the values 
of Pv/ Pc computed via the sums in equations ((47|) and 



(481 (open squares on Figure 17) to solve for pi/ Pc at low 



temperatures by "reflecting" them about the line deflned 
by eq. ([sT]). Thus 



(52) 



£i=2 + 2di.-^. 

Pc Pc 

The results are shown by empty squares on Figure [17] 
The error bars on pij Pc are equal to the error bars on 

Pvl Pc- 

B. Density 

To obtain a p\v-T coexistence curve in a non-reduced 
form (e.g. temperature in units of MeV and density in 
units of nucleons per cubic fermi) we flrst multiplied the 
temperature axis by T^. Errors on the temperature scale 
are then given by: 



8T = 6Tc 



(53) 



To determine the density in units of nucleons per cubic 
fermi we note that at T = the density of nuclear matter 
should be the density observed in unexcited nuclei. Using 
the same value of Tq = 1.2181 fm as in eq. (21 1, which is 



within 2% of other the leading order of other estimates 
[l8] , the density of nuclear matter at T = is 



PI {T = 0) 



Airrl 



0.132 A/fm^ 



(54) 



That value sets the scale on the density axis. The results 
are shown in Figure [TS] 

The error bars shown in Figure 



as the error bars shown in Figure 17 



18 



are not the same 
Simply translating 
the errors from upper right plot of Figure [iT] would give 
no estimate of the error of the critical density since, by 
deflnition there is no error associated with pd Pc- There- 
fore, we did the following: the error on a density value 
at a given temperature value p{T) is 



5p{T) 



\p{T + 6T)-p{T-ST)\ 



(55) 



Then we have a critical density value of pc = 0.06 ± 0.02 
A/fm'^ which agrees well with theoretical efforts [571 1571 - 
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FIG. 18: The temperature-density coexistence curve for bulk 
nuclear matter. Errors are shown for selected points to give 
an idea of the error on the entire coexistence curve. 



FIG. 19: The pressure-temperature coexistence curve for bulk 
nuclear matter. Errors are shown for selected points to give 
an idea of the error on the entire coexistence curve. 



Pressure 



To determine the coexistence curve for pressure as a 
function of temperature we again start with the reduced 
quantities and obtain p/pc as a function of T/Tc by per- 
forming the sums in equations ( 45 ) and ( 46 1 . We then 



determine the value of Pc from the compressibihty at the 
critical point which is defined as 



Zr = 



Pc 



PcTc 



(56) 



For fluids this is a universal quantity and is Zc — 0.277 ± 
0.004. Combing equations (46), (48) and (561 shows that 



the compressibility at the critical point is just a ratio of 
two Riemann ^-functions |50j 



Zr = 



C(r-l) 
C(r) 



=: 0.276. 



(57) 



However, using the error on r (given below Eq. ( 25 )) gives 
Zc = 0.28± 0.01. We use this value of Z^ in combination 
with the values of T^ and Pc determined above to obtain 
a value for the pressure at the critical point of 0.3 ± 
O.lMeV/fm'^ which agrees well with theoretical results 
[88t l89] . Here the error arise from the errors on Tc and 
Pc- Now to obtain the pressure in units of MeV/fm^ we 
multiply p/pc by the value of Pc obtained above. The 
error on the pressure is given by 



Sp = Spc 



P 



(58) 



Figure 19 shows these results. 



VI. SUMMARY 



The aim of this paper was to extract that liquid- vapor 
phase diagram of infinite, uncharged, symmetric nuclear 
matter from the data measured in various nuclear reac- 
tion experiments using finite, charged, asymmetric nu- 
clear matter, i. e. atomic nuclei. Because usual thermo- 
dynamical methods are obviously not accessible in this 
case, we concentrated on the fragment charge distribu- 
tions at various excitation energies and analyzed them 
according to Fisher's droplet model modified to account 
for the finite size of the system and the nuclear nature of 
the fluid (e.g. isospin and Coulomb effects). 

By fltting the charge yields observed in six different 
reactions studied in three different experiments, the crit- 
ical point was found to be: Tc = 17.9 ± 0.4 MeV, 
Pc = 0.06 ± 0.02 A/fm3 and Pc = 0.3 ± 0.1 MeV/fm^. 
Using the critical temperature and assuming the forma- 
tion of fragments exhausts all non-idealities, the entire 
coexistence curve of bulk nuclear matter was determined 
from T = to the critical point. This represents the first 
experimental measure of the phase diagram of bulk nu- 
clear matter and it is likely that the ideas and techniques 
outlined in this work would be useful in mapping other 
areas of the phase diagram of nuclear matter such as the 
phase transition between hadronic matter and the quark 
gluon plasma. 
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